Thermoelectricity by Perfectly Conducting Channels in Quantum Spin Hall Systems 



O 

U 

Or 
< 



i 

c 

o 
o 



(N 
> 

(N 

oo 
d 

On 
O 



X 



Ryuji Takahashi 1 and Shuichi Murakami 1 ' 2 

1 Department of Physics, Tokyo Institute of Technology, Ookayama, Meguro-ku, Tokyo 152-8551, Japan 
2 PRESTO, Japan Science and Technology Agency (JST), Kawaguchi, Saitama 332-0012, Japan 

(Dated: April 12, 2010) 

Thermoelectric transport of two-dimensional quantum spin Hall systems are theoretically studied 
in narrow ribbon geometry. We find that at high temperature electrons in the bulk states dominate. 
However, by lowering temperature, the "perfectly conducting" edge channels becomes dominant, 
and a bulk-to-edge crossover occurs. Correspondingly, by lowering temperature, the figure of merit 
first decreases and then will increase again due to edge-state-dominated thermoelectric transport. 

PACS numbers: 72.20.Pa, 73.43.-f, 73.50.Lw, 71.90.+q 



Thermoelectric conversion of heat into energy is one of 
the challenging topics in material science. The efficiency 
of thermoelectric energy converters depends on the trans- 
port coefficients of the constituent materials through the 
figure of merit. The figure of merit ZT is defined by 
ZT = aS K T , [I| where T is the temperature, a is the 
electrical conductivity, S is the Scebcck coefficient, and k 
is the thermal conductivity from electrons and phonons. 
Maximum efficiency of a thermoelectric conversion cy- 
cle depends on ZT, and the highest record of ZT is on 
the order of unity. It is an important but challenging 
issue to search for thermoelectric systems with larger 
ZT. There have been several proposals to overcome 
this conflict and to optimize the thermoelectric efficiency 
One of the proposals is the phonon glass and electron 
crystal 2](PGEC). Because the phonon carries heat but 
not charge, phonon conduction reduces thermoelectric ef- 
ficiency. Hence to achieve a high ZT, the system should 
be a bad conductor for phonons but a good conductor for 
electrons. These two conditions often conflict with each 
other, making materials search difficult. Another pro- 
posal is low-dimensionality. Q Low-dimensional systems 
have a peaked structure in the density of states, which is 
good for large S. Despite these proposals, good thermo- 
electrics have remained elusive and awaits qualitatively 
new approaches for improvement of ZT. 

In this Rapid Communication we propose that the 
quantum spin Hall (QSH) materials show enhanced ther- 
moelectric figure of merit at low temperature. The 
QSH systems are new state of matters for bulk insula- 
tors |J-|6|i realized in two-dimension(2D) and in three- 
dimension(3D). The 2D QSH system has gapless edge 
states which are stable against nonmagnetic impurities 
@, H|. Hence we expect that in dirty systems, elec- 
tron conduction through the edge states remain good, 
while phonon conduction is suppressed, satisfying the 
PGEC criterion. In addition, the edge states arc one- 
dimensional(lD), which fits the "low-dimensional" crite- 
rion. Another good reason for this expectation is that 
the QSH effect was observed in Bii^^Sb^Q, Bi9Se^[Io|. 
and Bi2Te3 ll| which are good thermoelectric materials. 

In 2D QSH systems in ribbon geometry, both the bulk 



states and edge states contribute. Because the number 
of bulk states is proportional to the ribbon width, we set 
the ribbon width to be very narrow, thereby the edge 
states can have comparable or even larger contribution, 
compared with the bulk. We then find that the bulk and 
edge contributions compete each other. We also find that 
there occurs a bulk-to-edge crossover when the tempera- 
ture is lowered. Because the edge states undergo inelastic 
scattering and lose their coherence, inelastic scattering 
length £i ne i gives an effective system size for quantum 
transport by edge states. As the temperature is lowered, 
And become longer, and the edge states become domi- 
nant in thermoelectric transport. We note that the edge 
transport cannot be dominant over the bulk transport at 
room temperature because £i nc i might become very short. 

The electric current j and thermal current w are cou- 
pled, and are induced by the thermal gradient or the 
electric field. In a linear response, they are described as 
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where q is the electron charge — e, and fj, is the chemical 
potential. Thermal and electric properties are given by 



a = e L , 
ZT 
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where n e is the electron thermal conductivitiy, and is 
phonon thermal conductivity. 

We first consider the edge transport only and neglect 
the bulk part. This corresponds to a case with very 
strong disorder, where the bulk states are assumed to 
be insulating, and the phonon heat transport is negligi- 
ble. To describe the coherent transport of the edge states, 
we use the Landauer formula. The density of states are 
schematically shown in Fig. [T] The edge states are as- 
sumed to be perfectly conducting over the whole sample, 
and the transmission coefficient T(E) is unity when the 
electron energy is within the bulk gap (—A < E < 0). 
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FIG. 1: (Color online) (a) Schematic bands for the bulk and 
edge states used in the calculation, (b) Thermoelectric figure 
of merit ZT as a function of chemical potential fi, by consid- 
ering only the edge states in the 2D QSH system. Transport 
by bulk carriers and phonons is ignored. 



Here we measure the energy from the bottom of the con- 
duction band, and A is the energy gap. To clarify an 
interplay between the bulk and the edge states, we focus 
on the bottom of the bulk conduction band and neglect 
the valence band. We restrict the chemical potential to 
be —A/2 <C \i. L v is given by 



(2) 



where the suffix e means the edge transport, h is the 
Planck constant. I and s are the length of the sample 
and the cross section of the sample. This is rewritten as 



sn 



(e* + I)' 



-dx 



(3) 



where p, = j^pf, and A = ■ We calculate ZT, by 
employing the gap size of Bi 2 Te3 (A =0.15 eV).The result 
(Fig. [I} shows that ZT becomes larger and well exceeds 
unity, when the chemical potential is in the bulk band. 
It is because the edge states carry large energy. 

In reality, when the chemical potential is in the bulk 
band, the bulk transport dominates, and reduces ZT 
from the otherwise large value. We treat the bulk and 
the edge transport independently, which is valid within 
the inelastic scattering length. We calculate the bulk 
transport by the Boltzmann equation as 



L h v = jdE (E-im)" (-|Q D(E) 



(4) 



where the suffix b means the bulk transport, and D(E) 
is the density of states, r is the relaxation time which 
is assumed to be constant. The bulk band is assumed 
to be parabolic with an effective mass m. For simplicity, 
we include only the first subband due to the confinement 
within the ribbon, by assuming that the gap between the 
first subband and second subband is large. The transport 
coefficients are then given by L v = L e v + L h v with 



b _ 4^2mk B Tn*c{k B Ty 

— : 



(e* + l) 2 



dx (5) 



where /x* is the mobility, and the coefficient c is the num- 
ber of the carrier pockets. 

We calculate these transport properties at T =1.8K. 
We again employ the parameters for Bi2Te3 as follows. 
The parameters for bulk transport is taken from those for 
bulk Bi2Te3. The electron effective mass is 0.02m e where 
m e is the electron mass, c is 6. (i* is measured at tem- 
peratures higher than 80K, and we estimated /i* to be 
2000cm 2 V~ 1 s~ 1 at T =1.8K by assuming that (i* satu- 
rates at lower temperature due to disorder. The effective 
system size £ is the inelastic scattering length £i nc \, and 
we assume I <~ l^tm, which is a lower bound of in 
HgTe quantum well at 1.8K [njj]. s is lOnm x 0.5 nm. kl 
is 0.1 Wnr'K -1 . which is expected from extrapolation 
from experimental data[I3| and theoretical estimate [l3j |. 
These parameters might have some error bars, because 
of the lack of the experimental data for Bi2Te3 thin film. 
The results are shown in Fig. [21 For these parameters 
the energy difference between the first and the second 
subbands is about 0.14 eV, and the chemical potential /j, 
is assumed be less than this energy. Many thermoelectric 
materials such as Bi2Te3 are narrow-gap semiconductors, 
and the effective mass is much smaller than the electron 
mass. Hence the subband structure is prominent, and 
the above assumption is satisfied without difficulty. 
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FIG. 2: (Color online) A calculation example of (a) Conduc- 
tivity, (b) Seebeck coefficient, (c) thermal conductivity, and 
(d) ZT as a function of the chemical potential. 

From Figs. [21 ZT has a maximum when the chemical 
potential fi is near the band edge. This results from 
a competition between the bulk and the edge states as 
follows. The Seebeck coefficient from the bulk states is 
larger when /1 is in the bulk gap, whereas that from the 
edge states is larger when /i is in the bulk band. Their 
effects tend to cancel each other, because their charges 
have opposite signs. Therefore, maximum of ZT occurs 
when fi is around the band edge. 
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For optimization of the thermoelectric figure of merit 
in QSH systems, we define the following dimensionless 
parameters from the prefactors in Eqs. (O and ([5]); 
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The parameter r represents the ratio between the edge 
and the bulk transport and g represents the ratio between 
the phonon heat transport and the bulk transport. These 
ratios characterize thermoelectric transport of 2D QSH 
systems. For each r and g we maximize ZT as a function 
of p. In Fig. [3j we show the maximum ZT max and the 
value of p = /l max giving the maximum. To focus on an 
interplay between bulk and edge transport, we restrict 
fi to be near the conduction band edge, and ignore the 
valence band, by putting A — > oo. From Fig. Ela), as 
a function of r, ZT max becomes minimum at r ~ 2.6, 
because of a competition between the edge- and bulk- 
state transport. This interplay is prominent in the plot 
of /Umax in Fig. [3jb) . The plot has a jump at around 
r ~ 2.6. As seen in Fig. 02c), at about r ~ 2.6, the plot 
of ZT as a function of p has two peaks, one from the 
bulk and the other from the edge. As r passes through 
2.6 from below, the peak from the edge dominates the 
peak from the bulk, and bulk-to-edge crossover occurs. 




FIG. 3: (Color online) (a) ZT max and (b) /i max as a function 
of r and g. (c) ZT as a function of p, for various values of r 
at g = 0.5. Bulk-to-edge crossover is seen by increasing r. 

In Fig. [3ja) we can see that at r = (no edge trans- 
port), the resulting ZT is sensitive to g, and it is impor- 
tant to reduce g by suppressing the phonon heat trans- 
port. However, disorder also suppresses electronic trans- 
port, and ZT is not enhanced so much. On the other 



hand, as r becomes larger, the result becomes insensitive 
to g. Disorder will enhance r, because the bulk mobil- 
ity becomes smaller. The ZT will then be enhanced. 
Generally, at low temperatures both r and g tend to in- 
crease as T decreases, as we explain in the following. As 
T is lowered, the mobility fj,* increases and eventually 
saturates, hl is given by kl — ^Cv^Il, where C is 
the phonon specific heat, vl is the phonon velocity, and 
li, is the phonon mean free path. As the temperature 
decreases, II becomes larger and saturates, while C de- 
creases; hence, kl first increases and then decreases at 
lower temperatures. From these behaviors, r and g tend 
to increase at low temperatures, possibly below around 
10K. An estimation using the above-mentioned param- 
eters for Bi2Te3 nanoribbon gives r = 9.4 and g = 8.2 
at T = 1.8K, which is located in the edge-dominated 
regime. We can estimate the crossover temperature for 
Bi2Te3 narrow ribbon taking into account the tempera- 
ture dependence fi* and kl in the similar way as in Fig. [21 
and assuming that t- me \, decreases as T~ 15 as has been 
observed in quantum Hall systems The crossover 

temperature is estimated to be around 5K-10K. 

To realize the edge-dominated transport, the ribbon 
width w should be much longer than the penetration 
depth A of the edge states, thereby we can ignore hy- 
bridization of the gapless edge states at the opposite 
edges. This hybridization induces a gap S ~ te~ w l x to 
the edge states 15], where t is the bandwidth (several 
eV). The penetration depth A depends on the systems, 
and in some systems such as Bi ultrathin film, it is esti- 
mated to be on the order of the lattice constant [l6| . As 
we set w = lOnm which is several decades of the lattice 
constant, the hybridization gap <5 is estimated to be on 
the order of mK. Thus in our temperature range above 
IK, this gap can be safely ignored. When we make the 
ribbon width to be much narrower, comparable to the 
penetration depth A, the edge states at opposite edges 
hybridize and opens a sizable gap [15j, killing the per- 
fectly conducting edge channels. 

In 2D QSH systems, elastic backscattering of edge 
states due to nonmagnetic impurities is prohibited 0, S] ■ 
Inelastic scattering is a key factor to characterize trans- 
port properties of the system. The electrons in edge 
states keep their coherence within the inelastic scattering 
length 4nei, which plays the role of the effective system 
size. We first estimate the electron-phonon (el-ph) in- 
elastic scattering length £i ne i, following the calculation 
on the quantum Hall (QH) system 17j. Here we as- 
sume the edge-state dispersion to be linear with velocity 
v c . We put the bulk wavefunctions to be proportional 
to sin(7ry/ui). By considering scattering by 2D longitu- 
dinal acoustic phonons, the relaxation time r is given by 
t -1 = (r ee )~ 1 + (r eb )~ 1 , where r ee , r eh are relaxation 
times by the edge-edge, the edge-bulk el-ph scattering. 
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Following Ref.[17| we obtain 

T ee l6pC 3 L V c ' T eb 



n 3 V 2 T 2 
pc\v c 



(8) 



ep is a screened el-ph scattering potential. If 

-19 T rp _ w „ = 10-6] i g/ m 2 ) 



where V e 

we take V ep = l(T iy J, T = IK, p 
cr, = 10 3 m/s, A = 10~ 10 m, W — 10~"m as an exam- 
ple, we get r ee - lCT 8 s, r eb ~ l(T 6 s. If v c ~ 10 6 m/s, 
^inci = v c t ~ 10~ 2 m. Experimental t- lue \ is much shorter, 
implying that el-ph scattering is not crucial among vari- 
ous inelastic scattering in the QSH system around IK. 

In addition to the el-ph interaction, the electron- 
electron (e-e) interaction also induces decoherence of edge 
states. There are two types of e-e interaction: edge-edge 
e-e and edge-bulk e-e interactions. The edge-edge e-e in- 
teraction is renormalized into the edge state action, and 
form the Luttinger liquid. Therefore this edge-edge e-e 
interaction does not cause dephasing if the system is clean 
enough and the edge channels remain perfectly conduct- 
ing well above Kondo temperature [18(. In disordered 
systems, it gives rise to a finite inelastic scattering time, 
while its estimate will be difficult. In addition, the edge- 
bulk e-e interaction also appears at finite temperature, 
and it depends crucially on the details of the system. 
Calculation of e-e interaction in the QSH systems is in- 
teresting but is beyond the scope of the present Rapid 
Communication. 

The inelastic scattering length ^ ino i is accessible experi- 
mentally. In the HgTe quantum well, nonlocal edge-state 
transport is observed [19( in 1/xm sample at 1.8K. It indi- 
cates that £j ne i is longer than the sample size, i- mc \ > Ipm 
at T = 1.8K. It is limited by the potential inhomogeneity 
due to gating. On the other hand, the inelastic scatter- 
ing len gth is measured in a QH system to be about lpm 
at IK [l4J, and is decreasing function of temperature. 
Based on these data we have used £- mc i — lpm at T — 
1.8K in obtaining Fig. [2 If the inelastic scattering length 
can be made longer, it will increase r and enhance ZT 
by edge-dominated thermoelectric transport. 

We address implications of our theory for 3D QSH sys- 
tems (topological insulators). Because the surface states 
on 3D QSH systems are not perfectly conducting, the ef- 
fect of surface states in 3D QSH systems on thermoelec- 
triciy will be less prominent than that of edge states in 
2D QSH systems studied in this Rapid Communication. 
Nevertheless, there can be one promising possibility also 
in the 3D QSH systems. In 3D QSH systems, protected 
ID states [20j of the crystal exist on line dislocations, 
depending on the bulk topological numbers. These ID 
states are perfectly conducting. Recently, a prominent 
magnetofingerprint was observed in a topological insula- 
tor Bi 2 Se3, and it is suggested that the phase coherence is 

It is also suggested 



retained over 2mm at around IK 21 



this scenario is true, they can be dominant in low tem- 
peratures, as we have shown in this Rapid Communica- 
tion. The estimated phase coherence length l- lnc \ ~2mm 
is three orders of magnitude larger than that we used 
in our calculation, and it is favorable for thermoelectric 
transport. 

Recently, an anomalous enhancement of the Seebeck 
coefficient at 7K is reported in p-Bi2Se3 22j. Though our 
2D model cannot describe three-dimensional p-Bi2Se3, 
we may attribute this enhancement to either surface 
states or ID states along line dislocations. In particular, 
the ID states form perfectly conducting channels, and 
will enhance the figure of merit. We note that in our cal- 
culation the edge and bulk contributions to the Seebeck 
coefficient has opposite signs, because the carrier charges 
have opposite signs (i.e. holes and electrons), and there- 
fore the Seebeck coefficient changes sign at the bulk-to- 
edge crossover by changing T. On the other hand, the 
Seebeck coefficient on p-Bi 2 Se3 does not change sign by 
lowering temperature. Within our interpretation this im- 
plies that the bulk carriers and the ID carriers have the 
same signs for the charge in the experiment. 

To summarize, we study thermoelectric properties of 
two-dimensional quantum spin Hall systems. The edge 
states become dominant in thermoelectric transport at 
low temperature, which might be below 5K-10K for nar- 
row ribbons. This bulk-to-edge crossover temperature is 
higher for longer inelastic scattering length of edge states. 
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21[, that the transport involved in this magnetofinger- 



print is carried by these ID states on dislocations. If 
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